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Abstract 

Using the zeta function regularization method we calculate the ground state 

energy of scalar massive field inside a spherical region in the space-time of a 

point-like global monopole. Two cases are investigated: (i) We calculate the 

Casimir energy inside a sphere of radius R and make an analytical analysis 

about it. We observe that this energy may be positive or negative depending 

on metric coefficient a and non-conformal coupling £. In the limit R — » oo 

we found a zero result, (ii) In the second model we surround the monopole 

by additional sphere of radius rg < R and consider scalar field confined in 

the region between these two spheres. In the latter, the ground state energy 

presents an additional contribution due to boundary at ro which is divergent 

for small radius. Additional comments about renormalization are considered. 
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I. INTRODUCTION 



Different types of topological objects may have been formed during Universe expansion, 
these include domain walls, cosmic strings and monopoles |J. These topological defects 
appear as a consequence of breakdown of local or global gauge symmetries of a system 
composed by self-coupling iso-scalar Higgs fields <3> a . Global monopoles are created due to 
phase transition when a global gauge symmetry is spontaneuosly broken and they may have 
been important for cosmology and astrophysics. The process of global monopole creation is 
accompanied by particles production 0. Grand Unified Theory predicts great number of 
these objects in the Universe fU but this problem may be avoided using inflationary models. 
From astrophysical point of view there is at most one global monopole in the local group of 
galaxies J3|. 

The space-time of a global monopole in a 0(3) broken symmetry model has been investi- 
gated by Barriola and Vilenkin ||. They have shown that far from the compact monopole's 
core the space-time is approximately described by spherical symmetric metric with an ad- 
ditional solid angle deficit. It is also possible to find solution for the Einstein equation 
coupled with an energy-momentum tensor associated with a pointlike global monopole. For 
simplicity we shall consider in this paper this singular configuration. In Ref. || a simplified 
model is presented in order to consider some internal structure for the global monopole. 

The analysis of quantum fields on the global monopole background have been considered 
in Refs. - @. It was shown, taking into account only dimensional and conformal consid- 
erations 0, that the vacuum expectation value of the energy-momentum tensor associated 
with a collection of conformal massless quantum field of arbitrary spin in this background 
has the following general structure 

where the quantities SI depend only on the solid angle deficit and spin of the fields. For 
scalar field this tensor was investigated more carefully by Mazzitelli and Lousto || and for 
massless spinor field by authors |§ in great details. 
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The above energy-momentum tensor has non-integrable singularity at origin therefore the 
ground state energy cannot be found by integrating the energy density. The same problem 
also appears for cosmic string space-time [ffDJ and in Minkowsky one with boundary condition 



on dihedral angle [II]. The calculation of ground state energy for the cosmic string space- 



time was considered in Refs. []r2| , |i"5| using different approaches. For infinitely thin cosmic 



string, specific global effect appears which leads to additional surface renormalization [|12 
The ground state energy of massive scalar field in the background of a cosmic string with 
internal nonsingular structure has been considered in Ref. [ 13| . It has been found that it is 
zero for arbitrary transverse diameter of the string. 

The nontrivial topological structure of space-time leads to a number of interesting effects 
which are not presented in a flat space. For example, there appear self-interacting forces on 



a massive point-like particles at rest. These forces have been investigated in Refs. ||T^JT5 
for cosmic string and global monopole space-times respectively. 

In the framework of zeta function regularization method (16| (see also [H]]) the ground 
state energy of scalar massive field can be obtained by 

E(s) = \m 2s U{s - i) , (1) 

which is expressed in terms of the zeta function £4 associated with the Laplace operator 
A = —A + ^TZ + m 2 defined in the three dimensional spatial section of the space-time. Here, 
the parameter M, with dimension of mass, has been introduced in order to give the correct 
dimension for the energy. In order to calculate the renormalized ground state energy we 
shall use the approach which was suggested and developed in Refs. fll8 H21|1 . 

In this paper we would like to discuss the ground state energy of scalar massive field in the 
background of point-like global monopole space-time inside a spherical region considering 
an arbitrary non-minimal coupling of this field with the geometry. Because the energy- 
momentum tensor has non-integrable singularity at the origin we would like to investigate 
two cases: (i) In the first we consider a point-like global monopole and calculate ground state 
energy using zeta function approach; (ii) in the second one, we consider a sphere surrounding 



the monopole and cut out internal part of it by an appropriate boundary condition for radial 
functions. This procedure permits us to reveal the role of the singularity In the limit of 
zero radius of inner sphere, this model corresponds to topological defects because there is 
no internal structure of monopole for arbitrary radius of sphere. 

The zeta function of Laplace operator on the point-like global monopole background 



has been considered by Bordag, Kirsten and Dowker in Ref. [p2 using the method given 



in Refs. []18|— pT|| . There, the general mathematical structure of zeta function and the heat 
kernel coefficients on the generalized cone have been obtained. Because the main emphasis 



of the present paper is on the ground state energy, we shall rederive in Sec.|TJ some specific 
formulas for our case which was not considered in Ref. [ p2fl . 

The organization of this paper is as follows. In Sec.0 we briefly review some geometrical 
properties about global monopole space-time which will be needed. In Sec.|TTT|, the zeta 
function of the Laplace operator on three dimensional section of a point-like global monopole 
space-time is developed. In Sec.jTV] we consider the zeta function for global monopole space- 
time cutting out by the sphere around the origin. In Sec.[V], the ground state energy of 
massive scalar field with arbitrary non-conformal coupling on global monopole background 
is considered for both above cases. In Sec.[VT], we discuss our results. The signature of the 



space-time, the sign of Riemann and Ricci tensors are the same as in Christensen paper |23 
We use units h — c — G — 1. 



II. THE GEOMETRY 

Global monopoles are heavy objects probably formed in the early Universe by the phase 
transition which occur in a system composed by a scalar self - coupling triplet field (f) a whose 
original global symmetry 0(3) is spontaneously broken to U(l). 

The simplest model which gives rise a global monopole is described by the Lagrangian 
density below 

L = \{d l( t ) a ){d l r)-\{rr-v 2 f ■ 
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Coupling this matter field with the Einstein equation, Barriola and Vilenkin || have shown 
that the effect produced by this object in the geometry can be approximately represented 
by a solid angle deficit in the (3 + 1) - dimensional space-time, whose line element is given 
by 

ds 2 = -dt 2 + a~ 2 dr 2 + r 2 (d6 2 + sin 2 6d<p 2 ) , (2) 

where the parameter a 2 = l — 8irr] 2 is smaller than unity and depend on the symmetry break- 
ing energy scale rj. The solid angle in the geommetry defined by (@) is 4na 2 , consequently 
smaller than Air. So this spacetime presents a solid angle defict given by Sfl = 32tc 2 7] 2 . We 
also can note that it is not flat. The nonzero components of Riemann and Ricci tensors, 
and scalar curvature are given below 



a 



2 



2(1 -a 2 ) 



For further application let us consider extrinsic curvature tensor on the sphere of radius R 
around the origin 

Kij = ViNj . 

Here Nj is outward unit normal vector with coordinates Nj = (0, a, 0,0). This tensor has 
two nonzero components 

e v R 



III. ZETA FUNCTION FOR POINT-LIKE GLOBAL MONOPOLE SPACE- TIME 



In order to calculate the ground state energy given by Eq.(|l]) we have to obtain the zeta 
function of the operator A in the neighborhood of point s = —1/2. For the calculation of 
zeta function we follow Refs JT^J20|22]. The zeta function of the operator A = —A + ^TZ+m 2 
is defined in terms of the sum over all eigenvalues of this operator by 
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in) 

Here A? n ) is the eigenvalue of operator B = A — m 2 . The eigenfunctions of the operator A 
defined in (^) which are regular at the origin have the form 

$( r ) = wA^ m (0^)j (A r ) , (3) 
V ctr a 

where Y\ m are the spherical harmonics and J M is the Bessel function of the first kind with 
index 



^=^ + ^) 2 + 2(l-^)(£-i). (4) 

A discret set of eigenvalues Xij can be found applying some boundary condition imposed on 
this function. Let us consider the Dirichlet boundary condition at the surface of a sphere of 
radius R concentric with the pointlike monopole 



%M-^R) = . (5) 
a 

Then, the zeta function reads 

1 oo oo 

c^-^) = EE(2/ + i)(AL + ^ 2 ) 1/2 - s , 

Z 1=0 j=0 

where the label R in the zeta function was introduced to indicate this kind of boundary 
condition. The solutions Ajj of equation (|) can not be found in closed form. For this 
reason we use the method suggested in Refs. |T^pO| which allows us to express the zeta 



function in terms of the eigenfunctions. According to this approach, the sum over j may be 
converted into contour integral in complex A-plane using the principal of argument, namely 

d , . „ _ , A 



(«(s - A) = J2(2l + 1) jf rfA 2 (A 2 + ™ 2 ) 1/2 - s ^ In A-"J M (£J2) 



where the contour 7 runs counterclockwise and must enclose all solutions of Eq.(0) on positive 
real axis. Shifiting the contour to the imaginary axis we obtain the following formula for 
the zeta function (see for details) 



C«(s - \ ) = £(2/ + 1) r dk(k 2 - m 2 ) 1/2 - s ^ In AT%( * fl) . (6) 

2 71 l=Q Jm OK a 

Here 1^ is the modified Bessel function. Let us use the uniform expansion for the Bessel 
function as below 



where t = 1/y/l + z 2 , r](z) = y/1 + z 2 + ln(z/(l + \/l + z 2 )) and 2; = kR/fia. The firsts 
coefficients and the recursion relations for higher ones are listed in ||24|| . This uniform 
expansion leads to power series over m, and the term un gives the contribution ~ l/m 3 ~ N . 
We shall make the calculations up to N = 3. In this case we obtain the following formula 
for uniform expansion of the logarithm of Bessel function 

ln(AT%(-i2)) = ~ z) - ~ ln(l + z 2 ) + -D^t) + ^D 2 (t) + ^D 3 (t) , (8) 

where 

D 1 (t) = Tx la t l+2a = -t~—t 3 , 
' ^ ' 8 24 ' 

D 2 (t) = V l2 a t 2+2a = —t 2 - -t 4 + -t 6 , (9) 
21 j z_, a,a 16 8 16 ' V ; 

w ^ ' 384 640 128 1152 

In the above expression we omit all constants which are not important for the calculation 
of zeta function. Adding and subtracting uniform expansion (Rl) in integrand of formula (^) 
we may represent the zeta function in the form 

1 7TI 2s jj 

<5(s - ~ 2 ) = N«(s) + (4T)V3r( ,_, ) E *(., *) , (10) 



where 



JV*( a ) = £(2Z + 1)^ r <k L 2 - ( A] 2 I (ii) 

x A jln/^/ix) - //77(x) + jln(l + x 2 ) - ^(t) - ^A^) - ^3°^)) > 



1/2-s 



A- 1 (s,R) = 52(21+1) 

a 1=0 

A (s, R) = -2n 3/2 mZ(0, s - -) 







A 1 (s,R) 
A 2 {s,R) 
A 3 (s,R) 



ft 



Z(0,s)-^-Z(2,s + l) 



7r 3//2 ma 2 



m * + ~) - 6-2(2, s + \) + \Z{A, s + ~) 



2(3 2 

irma 3 \ . 1062 , , 884 , 

" W r (0 ' s + 1} " "T" (2 ' s + 2) + T ( ' s + 3) 

1768 

-gg-W - + 4) 



(12) 

(13) 
(14) 

(15) 
(16) 



Here 2 -Pi = 2-Pi( — §' s — 1; |; — f^f )) is the hypergeometric function, /3 = mi? and 



Z(p,a) = T(q)52jl 



21 + 1 



cefj, 



The series in Eq.(|i~2"D 

00 

T(s) = 52(21 + 1) 



1=0 



1 (1 + «V//3 2 ) S V P 



r ( s - 1) P ocfi . 1, 

■2-Pl 7T r ( S - 7; 



(17) 



7T 



(18) 



can be expressed in terms of the same function given in Eq. flT7p. Indeed, one can use 
analytical continuation of the hypergeometrical function |24j 



, . 1 1 / Ofl 

2F1 1 ~,« - 



+ 



2' ' 2' 

r(l/2)r(l/2-a) 
r(-l/2)r(3/2 - a) 



/9 



1 + 



a/iT(l/2)r(s- 1/2) 

y r( s -i) 



1-s 



1, s — 1: s + 



V 



l+(f)' 



So, the first term in the rhs. of the above equation cancels the second one, divergent term in 
the sum (18) which is due to term in (||) (see Ref. [13fl). Now, one can use power series 
expansion for the hypergeometric function because its argument 1/(1 + (a/i//3) 2 ) is always 
smaller than unity, so we get 



K ' 2^ 1 1 'f^ T(n + s + 1/2) 



(19) 



Therefore in order to calculate of the zeta function we have to obtain an analytical contin- 
uation to the series Z(p, q). In fact we may consider only 
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Z ' -') = r ")§ (l + ovW - <20) 
because the other functions with p = 2,4,6,.. can be expressed in terms of 2(0, q) only. 
Substituting the value for \x given in Eq.(^) into Eq.(|20|) we obtain 



Z (°'^ 2F (^ 2 'S((7TWW' (21) 

where b 2 = (3 2 + 2(1 — a 2 )(£ — 1/8). This series is convergent for ^tq > 1. It is no difficult 
to obtain the analytical continuation of this series for small value of parameter b. Indeed, 
expanding 2 in powers of b we have 

oo (_-i\k 1 

2(0, s) = 2(3 2s Yl -T^ T ( k + s)b 2k ( H (2k + 2s - 1, -) . (22) 

For analytical continuation of this function in the domain 3lq < 1 and for great value of b, 
let us consider the series below 



oc 



This series, which has been considered in great detail by Elizalde [23], presents the following 
analytical continuation for great b 

FYs a b 2 ) ~ ^ V fclMl^W r 2/^4- V^s-1/2) 

— - — V n s ~ 1/2 cos(27rna)K s _ 1 / 2 (27rnfe) . 

r ( s ) 

Here is the Hurwitz zeta function and K n is the modified Bessel function. Differentiating 
this series with respect to a and putting a = 1/2 we obtain the analytical continuation that 
we need, which is the following 

~ 1 + 1/2 ^ ~ (-1)T(I+ S ) 

g (a + 1/2)2 + i , 2) .-^TT)+g i!r(s) 6 &(-l-2/,l/2). 
Taking into account this expression we obtain analytical continuation for function Z(0,q): 

z(o,8) |6 2 r( s -i) + 2 f til!r(/ + s )r 2 ^(-i-2/,i/2)| , (23) 



where 6 2 //3 2 = 1 + 2(1 — a 2 )(£ — 1/8)//? 2 . This function has simple poles for integer numbers 
q = 1,0,-1,-2,.... In order to obtain a renormalized value for the ground state energy 
we have to extract from our expression for zeta function flT0| ) the part which survives in 
the limit m — > oo. Moreover to calculate the zeta function up to degree m° we need only 
two terms from series ( |2"B| ) in which 1, 1/2) = 1/24 , £#(—3, 1/2) = —7/960 and three 
terms of T(s) which are given by Eq.(|i~9|). 

Putting this expression into Eq. (|i~9"D , and Eqs. (|l2|) - ([U|) and expanding over 1//3 = 
1/mR -C 1 and s, also collecting terms with similar degree on the mass m up to m° (we 
cannot here collect higher orders of m because we used uniform expansion up to this power) 
we get 



rn 



-2s 



+ 
+ 



( 47r )3/2 
7 

—iraR 



AtvR 3 
AttR 



m 



T(a-2) 



+ 



-2^' 2 R 2 



S(s 



m 



a 



(A 



12' 



; r( g -i; 
r(s-|; 



+ 



2vr 3 / 2 (A 



ira , . 1 



229 2 , 2tt . a2 1 7 . 

a 2 ) (A 2 --AH 

2520 ' aR y 6 240 ; 



Vis 



m 



Vis 



+ 



(24) 



Here A = 2(1 — a 2 )(£ — 1/8). All these terms are poles contributions in zeta function, all 
next terms will be finite for s — > 0. Comparing the above expression with that obtained by 
the Mellin transformation over trace of heat kernel (in three dimensions) 



1/2) 



1 



m 



, dtt s - 3/2 K(t) 

l(s-i)7o W (4tt) 3 / 2 



^m 4 



T(s 



T(s 



?■>'''',,,■"' " '''' 2 ' P^2 r(-S 1) R R pfi ^(s) , 



(25) 



we obtain the heat kernel coefficients : 
4vri? 3 



5, 



7 



3a 



5f = -2% 3/2 R 2 , Sf = -Trai? 

2 3 



( — 19) > 

a 12 



!i/9/ » In ttc^ / a 1 229 9 . 2tc , . 9 1 7 . 

2vr 3/2 A , B? = — A + a 2 A 2 - -A + 

v 12 ; ' 2 R v 12 2520 ; aiT 6 240 y 



(26) 



Those terms which are proportional to inverse degree of a come from exponential part of 
the uniform expansion given by (0), and respectively for T(s) (|T9"|). The terms which are 
linear in a 1 or a come from the series J2 u k/ H k in (Q). 
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Now we may compare our results with well-known formulas given in Refs. [26|,|22,[17],[27|]. 
The coefficients Bq, Bf, B^, B§ coincide with general formulas in three dimensions (all 

2 2 

geometrical quantities are given in Sec.|I|) 
AttR 3 



B 



R 



3a Jv 
B? = -2n 3/2 R 2 = 



dV 



(27) 



2 Jav 



dS 



B*= 7 -iraR-—(A-±-) = (l-A f KdV +\ { (trK)dS 
3 a 12 \o / Jv 3 Jav 



^ / (-96fft + 161Z + 8^ ifc iVW fc + 7(trK) 2 - 10 (trK) 2 ) dS . 

12 JdV 



192 Jav 

However some problems are connected with the term B2. The general structure of this term 
is the following (see P3| , P7| , for example) 



B* = [ b 2 dV + [ c 2 dS 
Jv Jav 



(28) 



where 



180 



is volume part, and 



1 n ik n lk + ^n mj n 



1 /1 



180 



1 /l 



(29) 



° 2 = I {I ~ e ) 11 {trK) + I {^6 ~ e ) ^ ~ ^0 UlkNlNk {trK) + ^TZ Ujk N l N k KV 



- —n u K il + — 

90 315 



I (trK) 3 - 11 (trK) (trK 2 '} + y (trK 3 



+ -n(trK), (30) 



is boundary contribution. Taking into account the results obtained in Sec.|| we have these 
terms in manifest form: 



b 2 = — j — to A 1 a 2 A A H 

r 4 4vr I V 12 120 / a V 6 240/ 



4a 2 



C2 



315i? 3 



(31) 
(32) 



We observe that the b 2 is proportional to 1/r 4 and the integral over volume in Eq. (2"8") 



will diverge at origin. This problem has already been discussed by Cheeger [281, Briining 
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and Seeley |29] and Bordag, Kirsten and Dowker [22] using partie finite of the integral. We 



regularize the expression for B 2 by restricting domain of radial integration 

f R dr f /. 1 17 „\ 2tt/ a9 1. 7\1 167ra 3 , . 

B 2 = - -to A + a 2 ) A 2 --A + H . 33 

Je r 2 \ V 12 120 J a V 6 240/ J 315R V ; 

After integration we take its finite remainder parts as e — > 0, and the expression obtained 

in this way coincides with that given in Eq. (p6|) . 

Our expressions for the heat kernel coefficients also agree with that ones obtained in 

Ref. f2"2fl . In that paper the heat kernel coefficients have been calculated for conformal 

case (£ = 1/8 in three dimensions). In order to compare both results we have to set 

£ = 1/8 (A = 0) in Eq.(|26|) and use the formulas of Appendix A from Ref. |22] for the three 

dimensional case d = 2. 



For renormalization which we shall discuss later we shall extract from zeta function (|i~0"D 



the asymptotic expansion (|24|) . Because all divergences at s — » are contained in fl24|) , we 
set s = in the remained part. After long calculation we arrive at the following formula for 
zeta function 



+ B ?™ 3 ^l + B ^m^ + B l m + b ?tBt> } • (34) 



2' \" 2' 



where 



B R (f3) = -Rm 4 B* - Rm 2 B R + RB R = -(3% - f3 2 b R + b R . (35) 

In order to exhibit the dependence on mass the m and on the radius of sphere R, we have 
introduced in above formula, the dimensionless heat kernel coefficients by relations 

b R = B R /R 3 = ^,b R = B R /R = 7 -na - —(A - i) , (36) 

b R = B R R = 7ra(A - — + —a 2 ) - —(A 2 - -A + — ) . (37) 
22 v 12 2520 J a K 6 240 ; V ; 

The function tt R (f3) tends to a constant for (3 and VL R (f3) = -B R ((3) In (3 2 + \^b R /2 /f3 + 

0(1/ j3 2 ) for (3 — > oo. The details of calculation and close form of fl R (/3) are outlined in 

Appendix |A|. 
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At this point we would like to make a comment. The origin of the term B R lnf3 2 is the 
following: In the limit m — > oo the singular part of zeta function has the structure given 
by Eq.(p5|). For small value of m it has the same poles structure multiplied by (3 2s . This 
is because all functions Z(p, s) are proportional to this degree of (3 as it may be seen from 
Eq. (P^D . The difference between them in the limit s — > is sin/? 2 multiplied by Eq.(^). 
Obviously that in this limit only B R , B r , B r survive which give the logarithm contribution 
to Eq.(|34]). 



IV. THE MODEL 

Because the geometrical characteristics of global monopole space-time are divergent at 
the origin we consider the following model: The center of monopole is surrounded by sphere 
with radius r whose interior region is cut out. It means that in our model there is no an 
internal structure for the global monopole. The present model reflects this peculiarity of 
topological defect. 

In frameworks of this model, we have to take into account both solutions of radial 
equation of the Laplace operator, instead of only one given in Eq.(^) which is regular at 
origin. The eigenfunctions now have the following form 

$(r) = J—Y lm {9^) \c 1 J lu {-r)+G 2 N^-r)\ , (38) 
V ar [a a J 

where is the Bessel function of the second kind. 

In this case we have two boundaries and one has to impose two boundary conditions. 
Let us again choose the Dirichlet boundary condition for the radial functions at spheres of 
radii R and r : 

C 1 J^-R)+C 2 N^(-R) = , (39) 
a at 

and 

C 1 J,(-r )+C 2 N,(-r )=0 . (40) 
a a 
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The set of discret eigenvalues Ajj can be found from equation below 

M^N^R) - N^r )J^R) = , (41) 
a a a a 

which is, in fact, the condition for existence of the solution Therefore we obtain the 

following formula for zeta function instead of Eq. (|6|) 

1 N COS ITS Cv 



a(a-i) = -— E(» + !) f ^ 

l IT l=Q Jm 



m 2 ) l/2- S 



9 , / T ,kR. Tr ,kr . Tr ,kR. T ,kr .\ . AnS 

x T^y^^ - K ^ )I ^ ) ) • (42) 

This general expression may be essentially simplified in the limit R/vq — > oo which we 
are interested in. Taking into account that in this limit the ratio K^ikRj <y)j l^ikRja) < 
7rexp(— ImRja) is exponentially small, so we may divide the expression for zeta function 
(|jD in two parts 

Us-h = c%(s-h+e A o (s-h, (43) 



2 y 2' 



where 

o 1 COS7TS 

Z 71 

and 



53(21 + 1) / dA^-m*) 1 / 2 -'^!* Ar%(-J2) , (44) 
Jm ok \ a J 



Ca°(s ~\) = -— + 1) f rfA; ^ 2 - "^"'/l ln f*^ V) • (45) 
2 7r /=Q Jm a/c \ a J 

The first part is the zeta function for pointlike global monopole which we have already 
calculated in last section. It depends only on the boundary condition on the sphere of 
radius R. The second part depends on boundary condition on the inner sphere of radius 
r . This kind of division of zeta function has been taken place for the case of thick cosmic 



string in Ref. |13[. It is also in qualitative agreement with [p^] . Indeed, according with ||20 
the internal solution gives Bessel function 7 M and the external solution gives function in 
expression for zeta function. The first part of zeta function (H3) depends on the solutions 
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which are internal with respect of sphere of radius R and the second part of zeta function 
( PD depends on the solutions which are external for sphere of radius r . 

Let us consider now the second expression (|45|). To calculate C A ° we use the same ap- 
proach which we have used in last section. We have to take into account the uniform 
expansion for modified Bessel function of second kind K^fxx) which has the form below 

W = (l + £(-1)^1 • (46) 

V r 1 I k=l P J 

Differently from the uniform expansion of Bessel function of first kind given by Eq. (]?]), the 
odd degrees of /i in above formula have the opposite sign. This fact leads to the change of 
sign of the heat kernel coefficients with integer index, also with respect of the heat kernel 
coefficients which were considered in last section. Using this uniform expansion we arrive at 
the following formulas for the zeta function 

1 m~'^ s _ 

Ca°(s --)= N^(s) + (4?r)3/2r(g - ^ J2(-l) k A k (s, ro ) , (47) 

where 

l/2-s 



v ; TlR 



Y(2l + l)na f°° dx | x 2 - ( — ) \ (48) 



x A hn(K^nx)) + m {x) + 1 ln(l + x 2 ) + -D x {t) - ^D 2 (t) + -^D 3 (t) 



and the functions v4fc(s,r ) are the same as in Eqs.(|T2|) - (|16|) but they depend now on the 
radius ro- Proceeding in the same way as it was done in last section, we obtain the following 
expression for second part of zeta function 

^ s - h = -isSa ln "° + °"<«} + { B ° m 'f(^T) (49) 

where /3o = ™~o and 

^ ro (/5o) = \r m 4 B r ° - r m 2 B[° + r B r 2 ° = ~ PX° + &2° ■ (50) 
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Contrary to the last section the heat kernel coefficients with integer number have changed 
the sign and they are 

4tt 7 4rc 1 

b r ° = BSVrg = -3^ , 6? = B?/r = --na + -(A - -) , (51) 

650 = B ?r = -vra(A - — + —a 2 ) + — (A 2 - -A+ — ) . (52) 
220 v 12 2520 J a K 6 240' V ; 

The summary heat kernel coefficients, according to Eq. (H3|) , are the sum of and B r n ° and 
they are in agreement with general formulas. We have to take into account that normal 
vectors for sphere of radius R and r have opposite directions and that the boundaries 
consist now of two spheres. It is easy to understand the division of zeta function in two 
parts given by Eq. (|43D , and opposite sign of the heat kernel coefficients B R and B r ° with 
integer indexes, by calculating B Q and Bi. For the space between two spheres we have 

* = I fl3 -I^< + B - < 53 > 

[ dS = -2n 3 / 2 R 2 - 27r 3 / 2 r 2 = B? + B? . (54) 

Jr 2 2 

Therefore the full zeta function in this case has the following form 

Us ~ \ ] = { Br ° {l3 ° ) ln0 ° + a '° {!3 ^ ~ { BR{I3) ln/32 + SiR(;3) ^ (55) 







4tt fR „ 


B = 


/ dV 


= — / r 2 


Jv 


a Jr 








Bi = 


[dS-^ 


2 


2 . 


Ir 2 



r(*-i; 



The close expression for Q ro is given in Appendix 0. 



V. THE GROUND STATE ENERGY 

In the framework of zeta function approach the ground state energy is proportional to the 
zeta function of Laplace operator and is given by Eq. ([]]). In order to analyze this energy, 
let us first of all, consider the ground state energy for point-like global monopole. The 
full energy of the system consists of two parts, namely classical part due to the boundary 
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and monopole background, and quantum one loop correction. The general expression for 
boundary contributions has been considered in Refs and it has the following form: 

E% = PrV r + a R S R + F R R + A R + ^-. (56) 

Here V R = 4vr J R 3 /3a and S R = AttR 2 are the volume and area of spherical surface, respec- 
tively. The two parameters p R and a R have simple physical means as pressure and tension 
of surface. The constant contribution described by parameter may be explained by the 
cosmological constant |3(J . The other two parameters F R , h R have not got a special names. 

The energy of monopole background can be obtained by integrating the (t, t) component 
of energy momentum tensor [/]] 

E 9 ™ = - ( R ^dV = -AixrfaR . (57) 
jo r 2 



The quantum correction, using Eq.(p4|), is 



E 



R 



2 bAV 2 ys u 32tt 2 /5 l r K ^ J i Km J 16tt 3 / 2 

where the heat kernel coefficients B R and B R are given by Eq. (|26|) and ([35"!), respectively. 

In order to obtain a well defined result for the full energy, we have to renormalize the 
parameters of classical part ( |56"D according to the rules below: 



' M \ 2s 3m 4 b R T (s - 2) m 3 b R /2 
Pfi^PR-l— IT ' a R ~^ ®R 



m , 



3m% R T 


> - 


2) 


64vr 5 / 2 r< 


s — 


1) 


m 2 6f T 


(*- 


1) 


16vr 3 / 2 T 




1) 








16tt 3 / 2 r 


— 


1) 



Fr ^ Fr -{^) i';l':V-> • a « a « - , (-™> 

h R ^h R - (—] 



96vr 5 / 2 ' 

mb i/2 
16tt 3 / 2 



V m / 

After this procedure we obtain the following expression for ground state energy 

Ef = ln/3 2 + . (60) 

The similar general structure of ground state energy in massless case has been obtained by 
Blau, Visser and Wipf [16] using the dimensional considerations only. For massive case we 
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find in manifest form the same structure. If we have used another scale for the mass like 
M — > M/x, in renormalization rules fl59"|), the above logarithmic term ln/5 2 will be replaced 
by \n(x(3) 2 . 

The expression ( |6"C]) is, in fact, the Casimir energy for internal part of the spherical bag 
in global monopole background. For small radius of the bag, this energy tends to infinity as 
In R/R: 

and for great radius of the bag R — » oo it tends to zero 

E « ~ "16^ • (62) 

Using these two limits we may analyze qualitatively the dependence of the Casimir energy 
of the internal part of the bag on its radius. The behavior of energy is defined by two heat 
kernel coefficients b 2 and bf, 2 . Both of these coefficients are the functions of non- minimal 
coupling parameter £ and metric coefficient a. In general, three kinds of different behaviors 
exist, which are plotted in Fig.l. It is possible to analyze the energy in general case, however 
we shall discuss only three cases for £ = 1/6, 1/8, 0. 

1. £ = |. In this case the behavior may be only of I and II kinds namely, the first kind 
for a < 1.24 and the second one for a > 1.24. The coefficient b 2 does not change its sign, 
but bf, 2 does for a = 1.24. 

2. £ = |. The behavior of energy may be I, II and III kinds, the first kind is for 
a < 1.045, the second kind is in the region 1.045 < a < 1.17 and the third one for a > 1.17. 
In the point a = 1.045 the coefficient bf, 2 change the sign, but b 2 does not up to a = 1.17 
where it changes the sign, too. 

3. £ = 0. This case is similar to previously one: it is of the first kind for a < 1.016, of 
the second for region 1.016 < a < 1.054, and of the third for a > 1.054. 

For a < 1 and £ = |, |, it is possible only the first kind of behavior. In the case when 



a = 1 the energy was calculated numerically in Ref. [20] and our results are in agreement 



with it. In this case bff = — 7r 2 16/315 and bf, 2 = — 7r 3 / 2 /120 the dependence may be the of 
first kind, only. 

In the limit R — > oo, the quantum correction tends to zero and the full energy contains 
only classical part which are due to boundary and background itself. 

Let us now proceed to our model. We surround the monopole origin by spheres of radii 
tq and R> tq and consider the bosonic matter field in the space between them. We do not 
take into account the interior of sphere of radius r , there is nothing inside it. We impose 
the Dirichlet boundary condition on this sphere which means that there is no flux into this 
region. The full energy in this case consists of five parts 

E = E* + E r J + E 9 ™ + Ef + Eg , (63) 

where 

E* = PbVr + *rSr + F R R + A R + ^ , (64) 

h r 

E d = Pr V ro + (T ro S ro + F ro r Q + A ro + — , (65) 

r o 

E 9 d m = -4nrfa{R - r ) , (66) 

are the classical part of energy due to the boundaries and global monopole itself, respectively, 
and 

K rn r Rfn ^ „ 2 , nR ,„^ , (M\ 2s 1 



E* = -M^(s - = [B*W In/? 2 + + % 



2 2 ,s ^' 32tt 2 /5 L ^' r V ^J \mj 16tt 3 / 2 

x { B ^W^f) - + ^W^T) + B l m + L ' (67) 

E r q ° = \m 2s Q{s - \) s ^ = {^(ft) ln/3 2 + Wfa)} + / }J ^ ' ' 



2 SAV 2 7S ^ U 327T 2 /? 1 1 1 \ m ) l 67 r 3 / 2 

x W^ T -^- - |i??m 3 + B^ T -^- + B?m + B?-£®\ , (68) 

I T( S -2) 32 ^-2) 2 r ( S - 2) J s ^0 

are the quantum corrections. Adopting the same renormalization prescription given in 
Eq.(|59~D for parameters in E^ and E^°, one arrives at the following expression for renor- 
malized quantum corrections 
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E R = -^L_ {B R ((3) ln/3 2 + n*(0)} (69) 

and 

^ ro = -g^r {^(ft) ln/3 2 + Q r °(Po)} ■ (70) 

The sum of these terms gives the Casimir energy of the field in space between the two 
spheres in the global monopole background. The first part we have already discussed. We 
may consider the second part in the same way. For small radius of sphere, r — > 0, it tends 
to infinity 

^~-!6^"' nA ' < 71) 
and for great radius of sphere r$ — > oo it tends to zero 

mb r ° 



-iro 



'5/2 



(72) 



Q 167r 3 / 2 /5 2 ' 

Because b^ = —b R and b r ^ 2 = b R 2 , the energy E r q ° has different behavior at small radius 
r . The sum of b R and b r £ constitutes the whole heat kernel coefficients for this space. For 
this reason the three kinds of dependences of E r ° on the radius are possible, only, which are 
displayed in Fig.l. The same results are available for E r q ° as it was obtained above for E R 
for £ = 1/6, 1/8, ; we have to change only the left plot to right one in Fig.l. The case 
62 = has to be analyzed numerically, however this discussion is out the scope of present 
paper. 

For any non-zero radius of the inner cavity r we have finite result. In this case the 
Casimir energy may be positive or negative, depends on the parameters of the theory. The 
main problem now is with the limit r$ — > 0, which has to reproduce the topological defect 
itself. The energy E r ° presents divergence in this limit as In r /r . This is in contradiction 
with above consideration of point-like monopole. If we set the radius r = at the beginning 



as it was done in Sec. Ill we obtain zero ground state energy for R — > 00. From the point 



of view of heat kernel coefficients we have already thrown away divergent part of B2 using 
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partie finite of integral ([33]). In frameworks of our model this thrown part appears here as 
divergent at origin and the additional renormalization is needed. 

For renormalization we may use the last term h ro /ro in the classical part of energy 
which is due to boundary (|5S|). We define a parameter Mq with dimension of mass by the 
relation h ro = GMq, where G is gravitational constant. With this definition, the divergent 
contribution for small radius ro may be canceled by the renormalization rule below 

2 

K ^K + ^h {2# ro ln(mro) + ^(mr )} , (73) 
where is the Plank mass, and then ground state energy is zero. 



VI. CONCLUSION 

In this paper we have considered the ground state energy of quantum scalar field in the 
background of global monopole space-time, with line element given by Eq.(0), in framework 
of zeta function approach. In order to reveal the role of singularity, we investigated two 
cases: In the first case we calculated ground state energy for point-like global monopole. We 
surrounded the origin by sphere of radius R and obtained that the ground state energy of 
the field inside. It has the form ( |60|) and tend to zero in the limit R — > oo and to infinity 
when R — > 0. The behavior of this Casimir energy in this cases is managed by two heat 
kernel coefficients and b§, 2 , respectively. The qualitative plots of the ground state energy 
for different values of the parameters £ = 1/6, 1/8, and a are given in Fig. 1(a). For a < 1 
it may be only in first kind for above values of £. 

In order to avoid the problem with singularity at origin, we investigated the second case 



in Sec.[IV], the following model: We surrounded the origin by a sphere of radius tq < R and 
considered the scalar field in domain between these two spheres using the Dirichlet boundary 
condition on the wave function associated with the massive scalar field, on the two surfaces. 
This boundary condition guarantee that there is no flux of particle through the spherical 
surfaces. The Casimir energy in this case consists of two parts given in Eqs.Q5"5|), Q7D|). The 
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first part is the same as for pint-like monopole case and second one is due to the inner 
sphere of radius r around origin. The structure of second part of ground state energy ( |70| ) 
is similar: there is logarithmic divergence at origin which tends to zero for infinite radius. 
The sign of energy for small distance is opposite. For this reason there appears three kind 
of dependence of energy which are displayed in Fig. 1(b). 

In the limit R — > oo and finite tq 7^ only one contribution in ground state energy 
survives (|69|). In the limit tq — > 0, it is divergent and additional renormalization is needed 
which is given by Eq.(fT^). After this renormalization the ground state energy of a global 
monopole will be zero. This is in agreement with ground state energy of point-like global 
monopole. 

If one fills up this cavity around the origin by matter, the situation becomes different. We 
may expect the same kind of divergence for additional energy from the interior of monopole 
but with opposite sign. We have already seen that the internal and external contributions 
have opposite signs: = —b^ - For this reason we may expect that this kind of divergence 
will cancel, however we cannot say anything analytically about divergence fl(mro)/rQ. In 
flat space-time |S{J it cancels, too, because the ground state energy is zero for zero radius of 
the bag. The same cancellation takes place in the case of a thick cosmic string background 
considered in Ref. fll3|. All of these aspects will be discussed in a separate paper. 
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APPENDIX A: 



In this appendix we want to give some brief explanation about the most important results 
found by us. First of all we represent the expressions for A k as series in powers of (5 — mr 
regarding for a moment b, (5 < 1 and A > —1/4 for convergence of series. They are 



a Asm- -b 2 °- t (- 1 )y r («- 1 +") v - 



A ( S , i2) = -^/J^tt 3 ' 2 f; { —^b 2n T(s + n- l)( H (2s + 2n - 2, 1 



?1=0 



TV. 



2" v 2 y ' 

T(s + | + a + n- i) 



(Al) 
(A2) 



/-* n=0 a=0 1 1 2 ttl J 



x g (i/j 2 + A) a 



1=0 "i 



2s+k+2a+2n-2 



(A3) 



In the above formulas we have used the following notations : v\ — I + \ , A = 2(1 — a 2 )(£ — 
1/8) and b 2 = (3 2 + A. As we can see only the first three terms in the expression for A-i , 
with n = 0, 1, 2, two terms in A , Ai , A 2 and one term in A 3 are divergent in the limit 
s — > 0. Extracting these terms we may set s = in the remaining series and we get the 
following result: 



m 



-2s 



(47T) 3 / 2 r(s 



-j-/? 25 {m 4 B*T(s - 2) + m 3 5f r(s - |) + m 2 5f r(s - 1) 
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Inl^pix) - pLrj(x) + - ln(l + x 2 ) - -D x - —D 2 - —D 3 



a 



7 7 1 / 11 

--C'(-3) - — + In 2 + A ( -2C'(-1) + ln2 

2 s v ; 160 240 V V ; 6 6 
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(A5) 
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Each of the above series may be analytically continued in terms of digamma function \l/ for 
arbitrary values of b and A. For example 
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This kind of representation is suitable for numerical calculations. 

Adding and subtracting the asymptotic expansion of zeta function (^) we obtain the 
following formula 
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2 J ' 16tt 2 /5 I ^ y yi > $ (4vr) 3 / 2 \ r(s - §) 

+ B ? mi W^) + B > 2 fj^|y + B f m + B "r^i) } • < A15 > 

where = Y^k=-i<^k{0) + ^f{(3). It is easy to see that the function fl B ({3) tends to a 

constant which may be calculated using the above formulas. Indeed, in the limit (3 — > oo we 
have to obtain asymptotic expansion of zeta function. Because one has already extracted 
the first five terms, so the next term will be -65/2- For this reason we get the following 

behavior for great (3 : VL B {(3) ~ — B R In (3 2 + \^Tb§//3 + ... . The coefficient 6f may be found 

2 2 

in the same way and it has the following form 

b § = = n^ 2 l-a 4 + - (A - -) a 2 - ( A 2 - -A + —)} . (A16) 

2 m 2 116 2 V 12/ V 6 240/ J V ; 

It is easy to obtain the formulas for zeta function from above expressions. The index k 

in Ak corresponds to a term which is proportional to [i~ k in uniform expansion of the Bessel 

function in Eq.(P). The uniform expansion of the modified Bessel function of second kind 



given in Eq.(46) may be obtained from uniform expansion of the modified Bessel function of 
first kind (^) by replacing the index fi to — \i. For this reason in order to obtain the formulas 
for the zeta function we may replace R — > r , /3 — > (3q and Ak(s,R) — > (— l) k Ak(s, r ) 
in the above formulas for zeta function £4. As the odd degrees of \i give contributions to 
heat kernel coefficients with integer index, they will change the sign. Therefore we have the 
following formula for the zeta function 



B r °m 

Y(s-l 
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+ 3l°™ 3 S^4y + 5[°m 2 ^4r + B?m + B?-^L\ , (A17) 
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where Q r °(f3 ) = EL-i(-l) W/?o) + ^°(A)) and 



u;?(/3) = 32ttV uJuf + A f°° , tfo, x 2 - - (A18) 



x ^ jln K^fix) + j27](x) + i ln(l + x 2 ) + -A - ^D 2 + ^D 3 
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Figure captions 

Figure 1. Three types of dependence of the ground state energy for the field inside the 
sphere of radius R (a) and outside the sphere of radius r (b). 
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